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The thermalization process is studied in an expanding parton gas using the Boltz- 
mann equation with two types of collision terms. In the relaxation time approximation 
we determine the criteria under which a time-dependent relaxation time leads to ther- 
malization of the partons. We calculate the entropy production due to collisions for 
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, the general time-dependent relaxation time. In a perturbative QCD approach on the 

other hand, we can estimate the parton collision time and its dependence on expansion 
time. The effective 'out of equilibrium' collision time differs from the standard trans- 
port relaxation time, r tr ~ (o? s ln(l/a s )T) _1 , by a weak time dependence. It is in both 

' cases Debye screening and Landau damping that regulate the singular forward scatter- 

ing processes. We find that the parton gas does thermalize eventually but only after 
having undergone a phase of free streaming and gradual equilibration where consider- 
able entropy is produced ("after-burning"). The final entropy and thus particle density 

, depends on the collision time as well as the initial conditions (a "memory effect"). Re- 

sults for entropy production are presented based upon various model estimates of early 
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Hard and semihard parton scatterings are expected to be the dominant processes in ultrarelativistic 
heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) 
energies J|,^] . These hard or semihard processes happen in a very short time scale after which a dense 
parton gas will be formed. However, this parton gas is not immediately in thermal and chemical equi- 
librium. Secondary interactions among the produced partons may lead to equilibrium if the interactions 
are sufficiently strong. Exactly how this parton gas equilibrate and what are the time scales are under 
intense investigations ®^ and that is also the focus of this paper. 

It is commonly assumed in relativistic heavy ion collisions that the matter expands hydrodynamically 
shortly after the nuclei collide and that it is in thermal equilibrium locally in space and time. In the 
classic Bjorken model JjJ| hydrodynamic expansion is assumed after a time tq ~ lfm/c. Yet, when initial 
parton formation times are estimated to be a fraction of a fm/c |3|,[5|,^| the expansion might be much more 
rapid than the typical collision time of the partons produced in the collisions. For times shorter than 
the typical collision time of partons immediately after the initial collisions, the rapid expansion is closer 
to free streaming than hydrodynamic expansion. Only at times much larger than the characteristic 
collision time may the parton gas thermalize, equilibrate and expand hydrodynamically 0. However, if 
the collision time increases with time the gas may never thermalize. The characteristic collision time is 
therefore a crucial parameter. Since it among other things depends on the density, which decreases in 
time, it may be time dependent. In fact, when the collision time is proportional to the expansion time, 



parton production. 

PACS numbers:12.38.Mh,25.75+r,12.38.Bx 
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the system neither expands hydrodynamically nor as free streaming at large times but somewhere in 
between fic|| . 

In this paper we study the thermalization process of a parton gas expanding in one dimension using 
the Boltzmann equation with two types of collision terms. In the relaxation time approximation we show 
generally for which time-dependent relaxation times the partons thermalize and estimate when they do 
so and how much entropy is produced in the collisions. Secondly, we describe the parton scatterings by 
perturbative QCD in order to estimate the important collision time and its dependence on expansion 
time. The effective "out of equilibrium" collision time differ from the standard transport relaxation 
time, Tt r ~ (a 2 s ln(l/a s )T) _1 , by a prefactor and a weak time dependence, which aer calculated in 
the Appendix, It is, however, still the Debye screening and Landau damping that screen the singular 
forward parton scattering processes. Finally, we will give results for entropy production based upon 
various model estimates of the preequilibrium condition. 



II. FREE STREAMING VS. HYDRODYNAMIC EXPANSION 



In order to have a tractable approach we assume that the spatial variations are sufficiently small 
along the longitudinal (z) direction so that we can describe the parton gas by the Boltzmann equation 



coll 



We also assume along the lines of the Bjorken model (?j that the transverse dimension of the parton 
gas formed in the nuclear collisions is sufficiently large that the initial expansion is one-dimensional. 
Furthermore, we assume the central rapidity regime is boost invariant in rapidity, i.e., depending only 
on the invariant time r = y/t 2 — z 2 . The Boltzmann equation thus reduces to Q 



/ df(p±,Pz,r) \ = ( df(p ± ,p z ,r) 
V dr J, \ dr 



(2) 



The collision term on the r.h.s. of Eqs. deter mines the equilibration and has been studied 

extensively within perturbative QCD [pj|^2|Jl^ -[l7[| near equilibrium in a quark-gluon plasma. In a 
relaxation time approximation, the collision term is 



f df(p±,p z ,r) \ = _ f - U 
{ 5r 



(3) 



where 9 is the relaxation time, which determines the time scale for thermalization and / eq is the 
equilibrium distribution function. 



A. The hydrodynamic limit 



When collisions are sufficient to thermalize the system, i.e., the corresponding relaxation times are 
short as compared to expansion times, 5<t, hydrodynamics applies and the distribution function in 
the local comoving frame is the thermal one 

/ oq (p) = (exp(^^)±l)- 1 . (4) 

If particle production is sufficiently rapid, chemical equilibrium can also be reached and fi ~ for 
gluons. Wc assume that the net baryon density is relatively low in the midrapidity region so that the 
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quark chemical potential vanishes as well. In the Bjorken flow model Ml the expansion is idealized in one 
dimension and assumed to be rapidity independent. Including viscous dissipation one finds generally 
in hydrodynamics jy] [l3| 

* + i±£ = in±f, (5) 

CLT T T Z 

where e is the energy density, P the pressure, 77 and £ the shear and bulk viscosities respectively. 
We assume that parton gas is a weakly interacting QGP consisting of relativistic quarks and gluons. 
Consequently, the gas is described by an ideal equation of state with pressure P = e/3 where the energy 
density is e = aT 4 ; a — ^§-(1 + %Nf) (a = 15.6 for Nj — /). A bag pressure can be included but will 
not change the entropy production that we are concerned with here. 

In the relaxation time approximation the shear viscosity is rj = 6e. In a weakly interacting QGP the 
temperature is the only scale besides factors of the interaction strengths, e.g., the Debye wavenumber 
is qo ~ gT. Thus the shear viscosity necessarily scales like 

V = VT 3 , (6) 

where fj is a dimensionless constant depending on the coupling constant g. The bulk viscosity £ vanishes 
for an ideal relativistic gas flifl . 

Solving Eq. (j^), the entropy density s = (4/3)e/T is obtained. The total entropy, S 7 is proportional 
to ts and is thus 0,|2) 

s = s »( 1 + ^ (1 - < ? ,2/S) ) 3 ' (7) 

where So, tq and To are the initial entropy, time and temperature. As seen in Fig. ^ the entropy 
production increases rapidly with viscosity fj. If the decoupling (freeze-out) or hadronization time is 
sufficiently late the entropy approaches the asymptotic value S = Sq(1 + fj / 2cit T ) 3 . 

To get an idea of the magnitude of 77 we refer to recent calculation within perturbative QCD Jl4|,[l(| . 
For weak couplings, a s ^,0.1, the quark and gluon viscosities depend on the coupling constant as 

fj~(a 2 Ml/a s ))-\ (8) 

where the logarithm arises from the sensitivity to screening of long range (small momentum transfer) 
interactions. The quark viscosity is approximately four times larger than the gluon because quarks 
interact more weakly. For strongly interacting plasmas, a s ^0.1, the Debye and dynamical screening 
must be replaced by some effective cutoff due to correlations in the plasma, which leads to fj cx a~ 2 . 
Within the range a s ~ 0.1 — 0.5, the quark (gluon) viscosity decrease from 77 ~ 50 (15) to fj ~ 1.5 (0.8). 
If the initial density of quarks and gluons, n^, is less than that in chemical equilibrium, n eq , then the 
viscosities are larger by a factor ~ n eq /rii. 
The corresponding viscous relaxation time 

r " = 5 7Tp- = Ta^ ' (9) 

is, however, larger than the expansion time and so hydrodynamics does not apply at early times as 
pointed out by Danielewicz and Gyulassy O]. The entropy production will be much lower and will 
be calculated in the following sections by solving the Boltzmann equation in the relaxation time ap- 
proximation. For comparison the upper limit on the entropy production in Eq.(|3^) from solving the 
Boltzmann equation with a finite collision time, 9 = t v , for the case 77 = 16 corresponding to 6 ~ 4ro is 
also shown in Fig. as a dashed line. Note that the entropy production continues long after and r n . 
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B. Free Streaming 



In the opposite extreme to the hydrodynamic limit examined above, collisions are absent and partons 
stream freely. According to Eq. ^ the distribution function evolves as || 

f(p)=Mp±,PzT/T ), (10) 

where /q is the initial distribution function at time To- 

As will be described in the following section the parton gas will stream freely until collisions thermalize 
the system around a time 9. The free streaming alters the distribution in phase space drastically. 
Initially the partons have large longitudinal momenta due to the high relative energy of the incoming 
nucleons in the nucleus-nucleus collisions and the relatively small transverse momentum transfer in 
hadronic collisions. However, in the one-dimensional free streaming expansion of the system at later 
times only those partons with similar longitudinal velocity will travel together locally in space and time. 
Thus the phase space separates the longitudinal momenta and the distribution function changes from 
a wide to a narrow one in p z locally in space and time (see Fig. ||). Collisions will then attempt to 
thermalize the system towards an isotropic distribution. 

When the longitudinal expansion has extended the system to a size similar to the transverse size, which 
happens at a time of order the nuclear transverse dimension r ~ R, three dimensional expansion takes 
over. Hereafter the densities will decrease rapidly reducing collisions drastically and a free streaming 
scenario is again likely. Around the same time, however, the system may break up, freeze-out and 
fragment. 



III. THERMALIZATON IN THE RELAXATION TIME APPROXIMATION 

Baym solved the Boltzmann equation in the relaxation time approximation Eq. (^|) with a relaxation 
time 9 independent of time || . He found that the parton gas started out free streaming and gradually 
thermalized to hydrodynamical flow on a time scale given by 9. The distribution function at any point 
in the local rest frame changes from being highly anisotropic, with only small longitudinal momenta, 
to being isotropic for r 9. A similar calculation by Gavin p"o[ ], however, with a collision time scaling 
linearly with time, 9 = at, gave a qualitatively different result. In this case the gas ends up somewhere 
between free streaming and thermal equilibrium depending on the values of a. For small a the collision 
time is short and the final state is close to Bjorken flow but for large a the state is closer to free 
streaming. The time dependence and magnitude of the relaxation time is thus essential for describing 
the degree of thcrmalization and its time scale. 

We will in the following solve the Boltzmann equation within the relaxation time approximation for a 
more general time dependence of the collision time proportional to the expansion time r to some power, 

9(XT P . (11) 

This covers the constant relaxation time of Baym (p = 0), the linear one of Gavin jl0| t c0 u = 1 jna ~ 
r (i.e. p = 1) as well as the near equilibrium transport relaxation time r tr ~ T _1 ~ n^ 1 / 3 ~ t 1 / 3 (i.e. 
p = 1/3). We shall study under which circumstances the parton gas thermalizes, estimate the relaxation 
times and predict the degree and time scale for equilibration. 

The solution to the Boltzmann equation (0) can be written in terms of an integral equation 
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/ = f (p^ Pz T/T )e- x + f dx'e x '- x f m (^p 2 ± + {p z r/T') 2 ,T',fi') , (12) 



where the time dependence of the temperature T and chemical potential /i in f eq (p, T, /i) are determined 
by demanding the energy density and number density (assuming no particle production) be the same 
for feq and / at any time, i.e., 

e(T, n) = J dT p E p f(p) = J dT p E p f cq (p, T, fj) , (13) 

n(T, l x) = J rfT p /(p) = J dT p f cq (p, T, fi) , (14) 

where dT p = d a p/(2ir) 3 . The function x(t) is related to r by 

, P^l, (15) 



x(t) = I dr'/e = 

J to r 



T T£ 

(r) 



where 6 = 6{tq). The marginal case p = 1 was studied in pli| . 

The evolution of the parton gas can be conveniently studied by taking moments of Eq. (|l2| ) with 
respect to particle energies. Summing over particle momenta gives a simple integral equation for the 
particle density which has the simple solution n(r) = tiqtq/t which also follows directly from Eq. ([s]). 
Multiplying by particle energy and summing over momentum we obtain 

e x g(r) = h{r /r) + f dx' e x ' g(r')h(T' /t) , (16) 
Jo 

where 

9(r) = ^% (17) 

and 



1 / sin" 1 v / T :r 
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h(r)= dcos(«)Vl + cos 2 (u)(r 2 - 1) = - [r + — . (18) 

Jo 2 \ Vl - / 

Here v is the polar angle of the particle momenta with respect to the z-axis. The function h{r) is a 
monotonically increasing function between h(0) = tt/4 and h(l) = 1. The function g(x) is calculated 
numerically and is shown in Fig. |^ for various 0$ and p < 1. 
Performing a partial integration on Eq. (hq) gives 



^'^[5(0^7^=0. (19) 

Note that by assuming an initially spherical symmetric distribution we differ slightly from Baym He 
assumes the distribution to be peaked in transverse directions initially. The term /i(ro/r) in Eq. ( |l6| ) is 
then replaced by unity and the r.h.s. of Eq. ( |l9|) changes to 1 — tt/4. The difference in initial conditions 
will not affect any of our subsequent arguments concerning the p-dependence of the thermalization. 

When p < 1, x(t) increases with increasing r as seen from Eq. (|l5|). The exponential factor e x thus 
weights large x' and Eq. ( |l9| ) implies for large x or t 

A-[g( T ')h(r'/r)} lT , =T = 0, (20) 

since rfr = 8dx. The slope of h(r) near r = 1 is 1/3 and so we conclude that g(r) oc r -1 / 3 for large r. 
Thus e(r) oc t~ 4 / 3 which is the one-dimensional hydrodynamical limit. 
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When p > 1, x(t) is negative and decreases with increasing r to a minimal but finite value x m i n = 
-t /<V(1 -p). By differentiating Eq. © we find that </(t) is always negative and therefore g(r) 
decreases monotonically. Yet, by definition, g is positive and must therefore have an asymptotic value 
larger than or equal to zero. From the integral equation (|l^) we see that this value at r = oo or 
x — a; max is non- vanishing and thus g(oo) > 0. Consequently, e(r) = g(oo)To/r at large times which is 
the one-dimensional free streaming limit. 

We conclude that thermalization will be reached when p < 1 and the parton gas will expand hydro- 
dynamically at large times whereas when p > 1 it will continue to stream freely. The marginal case 
p = 1 was studied by Gavin |l0[ who found that the parton gas ends up in a state between hydrody- 
namic expansion and free streaming depending on the size of the prefactor a where 9 = ar. When a is 
small the collision time is always relatively short and the gas equilibrates near the hydrodynamic limit. 
When a is large the collision time is always longer than the expansion time and the gas continues to 
stream freely. One should, however, keep in mind the finite decoupling, freeze-out, or hadronization 
time. When it is shorter than the collision time, which will be the case for large 8q or p close to unity 
or larger, the parton gas does not thermalize. 



IV. TIME DEPENDENCE OF THE RELAXATION TIME 



As we have just seen, the magnitude of the collision time as well as its dependence on expansion 
time is crucial for the equilibration. We shall therefore study the collision term and calculate it within 
perturbative QCD where we know the scattering matrix elements between quarks and gluons. We 
shall here only look at elastic scattering processes even though inelastic processes has been found to 
be as important in both the case of later parton chemical equilibration Q,|6| and energy loss of a 
fast parton going through a QGP [ [L9| . Weldon has, however, shown |2(| that inclusion of absorption 
processes as well as production and absorption of virtual particles reduces the energy loss further. Near 
chemical equilibrium emission and absorption should balance but the free streaming process drives the 
system away from equilibrium and thus inelastic processes is a possible candidate responsible for further 
equilibration. 

The collision integral for scattering particles from initial states 1 and 2 to final states 3 and 4 is 




X [/i/ 2 (1±/ 3 )(1±/4)-/ 3 /4(1±/i)(1±/2)]^(P1+P2-P3-P4), (21) 

where pi are the parton four-momentum. We assume they are massless, i.e., Ei = \pi\. The (1 ± /,) 
factors correspond physically to the Pauli blocking of final states, in the case of fermions, and to (induced 
or) stimulated emission, in the case of bosons, 1/2 is the statistical factor, 16 for gluons and 127V/ for 
quarks and antiquarks. |Aii2^34| 2 = \Mi2^34\ 2 / (16EiE2E 3 E4) is the matrix element squared summed 
over final states and averaged over initial states. For scattering of gluons 

quark-gluon and quark-quark interactions are just 4/9 and (4/9) 2 times weaker respectively near forward 
scattering (t ~ 0). In a t = lo 2 — q 2 channel, a singularity occures for small momentum q and energy w 
transfers. In a medium the t~ 2 singularity is screened as given by Dyson's equation in which a gluon 
self-energy Hl,t is added to the propagator 



G 



tr 



n 



L.T ■ 



(23) 



As was shown in jl4|-|L6| that Debye screening and dynamical screening due to Landau damping effec- 
tively screen the longitudinal and transverse interactions off in most transport problems at a length 
scale of order the Debye screening length q^ 1 . For small momentum transfer, q <C E\, E2, one can split 
the matrix element into longitudinal and transverse parts []l6| , 

2 



\M„ 



9 



1 



(1 — uj 2 /q 2 ) cos</) 



(24) 



_q 2 +Il L q 2 -uj 2 +U T 

where cos<^> = (vi x q) • (vi x q). The gluon self energies, Hl and Ht are given in the long wavelength 
limit ( q < T) by §1 



n L (<z,^) 

Il T (q, uj) 



1 



2q 



In 
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q — lu 



■(1 



i hi 



- to 



(25) 
(26) 



2q 2 ' 4q' q' 

The Debye screening wavenumber in thermal QCD is q 2 D = g 2 (l + Nf /6)T 2 where Nf is the number of 
quark flavors. 

The Boltzmann equation with the full collision term, Eq. (^l|), has been solved for quark-gluon 
plasmas near equilibrium and a number of transport coefficients have been calculated to leading orders 
in the coupling constant p4] , p^|j22| . For viscous and thermal relaxation as well as momentum stopping, 
the "transport relaxation time" is generally 



a>(-)AT 



(27) 



where A (the "fugacity" ) is the ratio of the actual density to the one in chemical equilibrium at tem- 
perature T . This relaxation time may be used at later times when the parton gas is near equilibrium. 
In Bjorken flow the temperature scales like T cx t" 1 / 3 , i.e., r tr cx r 1 / 3 . Since this power is less than 
unity, the parton gas should thus equilibrate according to the analysis in the previous section. 

In nuclear collisions the parton gas may be far from equilibrium when first produced and the expansion 
may also drive it out of equilibrium as in free streaming. Solving the Boltzmann equation thus becomes 
a very difficult non-linear problem that requires major computational efforts which is being undertaken 
in a number of parton cascade models as, e.g., in p). We take another approach in this paper. As 
mentioned above, hydrodynamics does not apply at early times because of long viscous relaxation times 



and the parton gas is expected to expand as free streaming initially. With this initial ansatz in Eq. (10) 



for the distribution function we can calculate the change in the distribution function at early times from 



the Boltzmann equation with the full collision term of Eq. (21). More specifically, from the entropy 
density 



s(t) 



p 



[/ln/T(l±/)ln(l±/)], 



we have calculated the entropy production 



(—) 

( dt )c 



E 



In 



1±/ 



(28) 



(29) 



in the free streaming phase by approximating / by /q [see Eq. (Al) in the Appendix]. The initial 



entropy production can be estimated analytically with the full collision term to leading logarithmic 
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order in the coupling constants and details of this calculation is given in the Appendix. Ignoring quarks 
the final result is [see (|A12|)1 



coll 



= — Ta 2 s T*\l g In 



ff 87T 4 



9r 



nXo, g a s To 



M 2 -) 

TO 



(30) 



where Ao, s = exp(/io,g/T)) is the ratio of the initial density to that in chemical equilibrium. 

We can match this entropy production to that obtained in the relaxation time approximation thus 
determining an efective and momentum averaged relaxation time 0. The relaxation time and entropy 
production during this early period, tq <C t <C 0, are different from later times t^9 when collisions 
change the free streaming distribution functions. The entropy production for the initial free streaming 
in the relaxation time approximation is from Eq. (|2£ 

relax 



coll 



fgb) 

OTo 
ne(T )/T 



fo - /cq\ \JV\ + (PzT/To) 1 



Mo 



To 



1 



rl ln(l + VI- r 2 /r 2 ) 



TO 
T 



TQ 



t>t . 



(31) 



As in Eq. (J3^) we ignore quarks; e(T ) 
equating the entropy production of Eq. 
term Eq. © wc find 



a ^ g (7T 2 /30)To A ,g is the initial gluon energy density. By 
5fi) for r ^ tq to that from the perturbative QCD collision 



- ~ 1.43a 2 T Ao,(, In 



37TT 



2X , g a s To 



K— ) • 

TO 



(32) 



Note that it is the initial temperature To that enters here and not T(t) as in Eq. (|27|). The collision 
time depends only logarithmically on the expansion time. As explained in the Appendix, the relaxation 
time is only weakly time dependent in a free-streaming parton gas because the phase space for small 
momentum scattering opens up quadratically with time thus effectively compensating the decrease 
in parton densities. On the other hand, if large momentum transfers are imposed to each parton 
scattering, the entropy production rate will decrease quadratically with time, leading to a much stronger 
time dependence of the relaxation time. The long range interactions (small momentum transfers) are 
therefore very important in an expanding parton gas. The two logarithms in Eq. (^2|) arise from 
integrals over momentum and energy transfers respectively. The "fugacity" factor Ao arises from the 
correspondingly smaller density of scatterers. 

To estimate the relaxation time 9, we need to know the initial condition To, To, Ao = e M °/ T ° which so 
far can only be determined from model calculations. From the HIJING model calculation |,||J2| it 
was found that at r = 0.7 fm/c the produced partons in central Au + Au collision at RHIC energy can 
reach a local isotropy in momentum distribution temporarily with effective temperature To = 0.57 GeV; 
in addition X . g = 0.09 and A , g = 0.02. At LHC energies r = 0.5 fm/c, T = 0.83 GeV, A , 9 = 0.14 
and Ao.q = 0.03. Eskola et al. [g5| found similar results for the temperature To = 1(1.5) GeV at the 
RHIC (LHC) energy in their minijet plasma calculation. However, they used a smaller initial time, 
To = 0.1 fm/c, and find Ao l9 ~ 1 and Ao, g ~ If one were to allow this minijet plasma to stream freely 
shortly after To the result would be consistent with HIJING estimates at the later time To ~ 0.7 fm/c. 
The newest set of parton distribution functions in the calculation of Eskola et al. also increase the 
initial parton density, especially at LHC. In the parton cascade model Q the initial parton density is 
found to be larger due to a different treatment of soft parton interactions. These numbers for tq and 



8 



To are surprisingly similar to those found by Shuryak || in a different analysis. By estimating the 
particle rapidity distributions dN/dy in relativistic nuclear collisions he obtains a particle density in 
the Bjorken scenario (dN/ dy) / (ttR 2 t) which at a time r = r co ii ~ (a s T) _1 is assumed to be the same 
as the equilibrium one ~ T 3 . Hereby the initial values for T and r are found. One should bear in mind 
that all these estimates are based on perturbative QCD inspired models. There are many uncertainties 
(see, e.g., p6|) due to our limited knowledge of strong interactions. 

If we ignore the slow logarithmic time dependence and assume To — l/ T o an d ^o,q ~ Oj as motivated 
by the above mentioned models, we find 

9 ~ c*7 2 Ao> • (33) 
With a s ~ 0.3 we find a rather long collision time even if Ao, s = 1. 



V. ENTROPY PRODUCTION 



With the weak time dependence of the relaxation time, we can now estimate the entropy production 
during the early thermalization. We still assume the relaxation time has a power dependence on time 
as in Eq. (Ill]) with small p. The scaling behavior of the function g(r) with 9q — 6(t = tq) can be 



obtained when #0 3> Tq. In that case the integral equation (16) for g(r(x)) depends almost solely on x, 
since x' /x ~ (t'/t) 1_p . Thus for given p, g(r(x)) is a generic function of x. Its behavior at large x is 



g(r(x)) ~ gpar-Vsfi-P) 



fJ P 



TO 



1/3(1-P) / T \ -1/3 



TO 



(34) 



where g p is some p-dependent constant of order unity. In the case p — Baym [H found g p = 1.22 but 
due to the different initial conditions (isotropic versus peaked in the transverse directions) this value 
differ from our case by a factor 7r/4 when 9o ^> To. 

Until the parton gas reaches equilibrium its entropy at time r is always less than the equilibrium 
entropy at temperature T(t), i.e., s < s eq — (4/3)e/T = so5(t)ToTo/(Tt). Thus, the total entropy is 

S<S g(r)^=So (£) V4 girfl* . (35) 
If the parton gas equilibrates the equal sign holds at large times and the final entropy is from Eqs. 

(EE® 

^ = * p (Vto) 1/4(1 - p) , (36) 

where s p is a p-dependent number. In Fig. ^, the final entropy is plotted for various values of p and 0q. 
The formula Eq. ( |36|) is a good approximation with coefficient s p ~ 1. For comparison the final entropy 
calculated by Baym 

^ = 1.16(Vro) 1/4 . (37) 

with a constant relaxation time (9 — 9q, p — 0) is also shown in Fig. |], being multiplied by a factor 
(7r/4) 3 / 4 because of the different initial conditions. 

It is interesting to compare the entropy production to the case when the parton gas is streaming freely 
until a time r su dden when the gas suddenly collides violently and immediately thermalizes. Conserving 
energy per volume (but not particle density) we find 
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^ = (j) 3/4 (r S uddc„/r ) 1 / 4 , (38) 

which is only slightly lower than Eq. ([37]) if 9o — T su dden- The exact way by which the parton gas 
equilibrates is therefore not so important; it is the collision time that determines the entropy production. 

From Eq. (|37| , ^6|) it is evident that the final entropy production increases slowly with the collision time. 
Yet it takes much longer time and at a given time (t?S9) the entropy production rate [cf. Eq. ([Si])] is 
inversely proportional to 9. If the parton gas decouples, freeze out or hadronize at time, Td, and entropy 
is no longer produced, then the total entropy production will decrease as 9 increases above r^. 

With the approximate collision time of Eq. (|3^) and assuming To ~ 1 /to the final entropy is obtained 
from Eq. (§6|) or @ 

| -«; 1/2 AoT- (39) 

For a s ~ 0.3 this gives an increase in entropy by a factor of 2-3 when varying Ao from unity down to 
0.1. 



VI. SUMMARY 



We have studied a one-dimensionally expanding parton gas created in the wake of nuclear collisions. 
Within the Boltzmann equation in the relaxation time approximation we find that the rapid expansion 
is closer to free streaming than hydrodynamic expansion for times shorter than typical collision times 
of partons, t <t<9. Only at times much larger than the characteristic collision times, r 3> 9 may the 
parton gas thermalize and expand hydrodynamically. However, if the collision time increases with time 
the gas may never thermalize. Parametrizing the collision time as 9 = #o( r / r o) p the condition for 
equilibration and hydrodynamical expansion is p < 1. We calculate how much entropy is produced in 
the collisions. These different expansion modes and the corresponding parton phase space distributions 
at different stages of evolution in high-energy heavy ion collisions are illustrated in Fig. |[ 

We have described the parton scatterings by perturbative QCD in order to estimate the important 
collision time and its dependence on expansion time. The effective "out of equilibrium" collision time 
differ from the standard transport relaxation time, l/r tr ~ cq \n(l/a s )T(r), with a weaker time depen- 
dence. However, it is still the Debye screening and Landau damping that screen the singular forward 
scattering processes. 

We find that the parton gas does equilibrate eventually with these collision times but only after 
having undergone a phase of free streaming and gradual thcrmalization where considerable entropy is 
produced ("after-burning"). The final entropy and thus particle density depends on the collision time 
as well as the initial conditions (a "memory effect"). For various models predicting the preequilibrium 
scenarios the entropy production is significant. The total entropy and particle production is estimated 
to be doubled or tripled with respect to the initial value. 

These estimates do not include particle production which by itself adds to the entropy production. 
On the other hand particle production will also increase the density and thus shorten the effective 
collision time which leads to a decrease in entropy production according to Eq. (^) . 

Most analyses assume a constant density in space but large density fluctuations may well be present 
in the initial parton plasma. This will increase the average entropy production for both elastic and 
inelastic scatterings since these are proportional to the initial densities squared as well as the final 
densities through the stimulated emission factors (for bosons) or Pauli blocking factors (for fermions). 
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High density regions ("hot spots") will equilibrate thermally and chemically faster than low density 
regions. At the same time, however, the free streaming will tend to reduce density fluctuations. 

It was emphasized in the Appendix that the very singular small momentum transfers provides strong 
scattering and the opening up of phase space compensates for the decreasing densities. If a larger 
momentum transfer cut-off of the order of particle momenta (~ temperature) is applied then the collision 
time will increase quadratically with expansion time and the parton gas will never thermalize. Also, 
when the longitudinal extension of the system exceeds the transverse size, which is in the order of the 
nuclear radius, the expansion proceeds in three dimensions. The densities will then decrease cubically 
with expansion time and collision might never catch up with the expansion. 
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APPENDIX: THE RELAXATION TIME 



In this appendix we give a detailed derivation of the entropy production by elastic parton collisions 
in a system near the free streaming case at times tq <C t <C 9. 

The initial distribution of partons has been estimated in several models |^,^3|]. Typically, one finds 
that the partons are formed within the first tq ~ 0.2 — 0.7 fm/c after the nuclear collision and that rapid 
longitudinal expansion takes place. The local momentum distribution of partons at the formation time 
is not isotropic but forward/backward peaked, i.e., (\p z \) S> However, due to streaming the local 

distribution changes rapidly since (\p z \) scales as tq/t [see Eq. (|l0|)] and at later times the particles 
have (\p z \) <C (p_i_). At the cross over time To the parton gas is isotropic in momentum space and seems 
to be in approximate thermal equilibrium with temperature To even though it is streaming freely in 
space and time. From Eq. (|l0|) we obtain the free streaming distribution function 



/o(p,r) 



exp(-^^ — ) ± 1 



(Al) 



where /Lto is the chemical potential determined by the density. Since the parton gas may not have reached 
chemical equilibrium yet, /xq may not vanish. In fact most models predict that particle densities are 
rather low initially j|,|24|,|25|] and that — /xq — (1 — 2)Tq. The temperature increases with collision energy 



and typically To — 0.5 — 2GeV. We shall use Eq. ( |10[ ) with Eq. (Al) for the free streaming initially 
with parameters To, To and /io. The entropy production due to collisions is negligible at times around 
To because the parton gas is near thermal equilibrium and so we shall ignore collisions earlier than tq. 
On the other hand, continuing particle production will produce entropy but we shall not consider that 
contribution here. 

With the free streaming distribution function, we obtain from the Boltzmann equation by changing 
variables from p z ,i — > p z ,iTo/T, i=l,2, and q z — > q z T /r 
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'^1 

dt. 



-2-kv x vi T -^ J dr pl dr P2 dr q |Mi 2 ^34(q±, q z ^)\ 2 fo(E 1 )f (E 2 )(i ± / (£ 3 ))(i ± h(E 4 )) 

E\ - 



[1 - cxp( 



£1 



)}5(E 1 + E 2 -E 3 - E A ) , 



(A2) 



where £3 = y/ (pi± + q±) 2 + {pi, z + Qz) 2 Tq/t 2 and E\ = yp[j_ + {p z to/t) 2 . The same expressions are 
valid for i—2 : 4 when the sign of q is changed. 

At this point we want to emphasize the importance of small momentum transfer processes, q ~ 
Id ~ gT as compared to large momentum transfer ones, q ~ T. For the latter the exponential in 
Eq. (A2) can be ignored because £4 + E 3 — E 2 — E\ ~ T. One then finds that the integrations over 
the energy conservation 5- function removes a factor To jr leaving an entropy production rate decreasing 
quadratically in time. This just reflects that the particle densities of the two scatterers decrease as tq/t. 
The small momentum transfers have, however, a very singular scattering matrix element and, as we 
will now show, the phase space opens up quadratically with time for small q - effectively compensating 
the decreasing densities of scatterers. 

Expanding around small q to second order, the term in the square bracket representing the difference 
between scattering in and out is 



1 — exp( 



Ei + E 2 - E 3 - £4 (vi - v 2 ) • q [(vi - v 2 ) • q] 2 



q 2 - (vi • q) 2 q 2 - (v 2 • q) 2 



2T 2 



(A3) 



2-EiTo 2E 2 T 

Only when r = tq is E = E and energy conservation requires that Eq. ( |A3j ) vanishes. Due to symmetry 
the first term vanishes when integrated over pi and p 2 and the leading term is second order in q 2 . 
For small q we can also replace £3 by E\ and £4 by £ 2 in the distribution functions. Thus we have 



= — 2-kv\v 2 - 



coll 



dr p j (E 2 )(i±fo(E 2 )) / dr 



Ei - Mo 



pi 



T, 







/o(£i)(l±/o(£i)) 



dr q |Mi 2 ^34(qj_, <lzj)\ 2 j d ^{u -Ek + E 3 )S{uj + E 2 - £ 4 ) 
q 2 - (vi • q) 2 q 2 - (v 2 • q) 2 [(vi - v 2 ) • q] : 



2£!T 2£ 2 T 2T 2 

Here we have introduced an auxiliary integral over energy transfer u>. 

We can use up these S- functions by performing the angular integrals dQi, i 

,- 2 TT 



(A4) 



=1,2. For example, 



h 



dn 1 s(uj-E 1 + £ 3 ) = 



'1 ; sin^id^i 
'0 



x 5 



cos 0i sin 0%q± + cos 8iq z TQ /t 2 



(A5) 



V1-COS2 ^(l-T 2 ^) 

where 9\ and (f>% are the polar and azimuthal angles of pi in a coordinate system with z-axis along the 



collision beam direction. For t — tq the prefactor in Eq. (A3) vanishes due to energy conservation. For 
t > ro we can ignore the (<7zTo/t) 2 term in Eq. (A5), which then yields 




y 



q±_ Vi - bx 2 



(A6) 



where b 



tZ/t 2 , 



= cos 0\ and y 
— cj 2 /g 2 )/(l — buj 2 /q 2 J, this integral gives 



cos 4>i- Changing variables to sinx 
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h = r. \ n i , / 



tt/2 



X 



2 



4 , for to/t « 1 . (A7) 

qx\/l~buj 2 /qj_ 



The angular integrals of dD,2 yields the same factor as Eq. (A7) and the integral over energy transfers 



in Eq. ( |A4| ) thus gives at large times 

J duhh = J duidQ^QiSiui -E x + E 3 )5(u + E 2 - E A ) 



du 32 , / It 



16/ 7-^-— In — • (A8) 

J -ii. Tx - hbJ U \ T o. 



An angular dependence arising from Eq. (A3) should also be included when performing integration 
in Eq. (A5). For example, an extra factor of cos 2 0\ leads to an additional factor of 1/3 in Eq. (|Acf). In 
addition the matrix element depends on energy transfer through the transverse part of the self-energy 
in Eq. (^4j). After evaluating the integrals over momentum transfers this dependence is, however, only 
logarithmic in oj Jl6| and can be ignored. 

Let us first consider the entropy production due to gluon-gluon scatterings and include quarks later. 
We assume for convenience that ^to^ — Tq which allows us to use Boltzmann distribution functions, 



fo(E) = exp((/io — E)/To) = Ao exp(— E/T). The momentum integrals of pi and P2 in Eq. (A4) are 
straight forward, which leave one remaining integral over momentum transfer, 

'ds\ 3v 2 n ,r 3 , ,2r, 3 f .„ q 2 



at 



coll 



where, we have for simplicity approximated 11^ by q 2 D in the matrix element of Eq. (|22]). Furthermore, 
we have replaced the transverse part of interaction by the longitudinal one times a factor (cos 2 (f>) = 1/2. 
Both these approximations are exact to leading logarithmic order in the coupling constant for the 
calculation of a number of transport coefficient Note that the limit for the integral over q z is 

iftnaxT/To and the maximum momentum transfer, g max , is determined by the distribution functions 
which cut off large momentum transfers. It has been estimated in [fl6[ to be q max ~ 3Tq. The integral 
over momentum transfers in ( |A^ ) then gives (t/to) 3 7t 2 \n(2q max / qo) when r/ro (jmax/lD- It includes 
the usual logarithm (see, e.g., [|l6| for details) as well as a factor (t/tq) 3 from the integral over q z . The 
entropy production in QCD to leading order in the coupling coupling constant is thus 

2\2 ,4 \ \ i f^\2 2\2 i„ 4<? n 

v g A o,g + -Q V g v qAQ,gAQ,<i + \g) v q \), q 



ds\ 9 



4„2 



dt J coll 8^ 4 



In 



imax 



Q 2 n 



ln(-). (A10) 

TO 



Here, we have included contributions from quarks and antiquarks to the entropy production; v g = 16, 
v q = 12Nf, \ q = e^ q / T and we assume \ q — \q. The quark-gluon and quark-quark forward scattering 
interactions are smaller than the gluon-gluon ones by a factor (4/9) and (4/9) 2 respectively. According 
to the models in Refs. fewer quark and antiquark than gluons are produced in relativistic nuclear 

collisions, i.e., \o t qq *C Ao. g . 

The Debye screening mass in a quark-gluon gas in thermal and chemical equilibrium with no net 
baryon density (i.e., fj, g = (i q = fig = 0) is qfj — 47r(l + Nf/6)a s T 2 . Out of chemical equilibrium 
when [i q — [iq and /i g . g <^ —T the Bose and Fermi distribution functions can be replaced by Maxwell- 
Boltzmann distribution functions and we find 

q 2 D = -(\ 9 + ^fX q )a s T 2 . (All) 
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The quark and gluon densities and thus fugacities decrease with time due to the one-dimensional 
expansion such that A = XqTq/t. 



Inserting these expressions for g max and qo in Eq. (A10) we find the entropy production rate 



coll 



9 ™2 T 4 



2\2 i ^ \ \ / * \2 2 \2 



x In 



3ttt 



2(A , 9 + ¥fXo, q )a s To 



H 2 -). 

TO 



(A12) 



to leading order in a free streaming gluon gas. The result is approximately valid for t^2tq and as long 
as the free streaming assumption is valid, i.e., t^8. 

We emphasize the important result that r 3 factor is cancelled by the integral over momentum trans- 
fers. Only a slow logarithmic dependence on time remains. The physical explanation for this cancellation 



is the following. From Eq. ( A2) a factor (to/t) 2 appears from the substitutions of integration variables 
p z ,i,i =1)2. This represents the fact that the densities of each of the colliding partons drop like t /t. 
If we keep the original momentum transfer variables (i.e., do not replace q z — > QzTo/t in Eq. (A2)) 
the factor in the square bracket, Eq. (|A3|), leads to a factor (<7t/t ) 2 . This is because the phase space 
for small momentum scattering opens up quadratically with time and it balances the decrease in partem 
densities. On the other hand, for large momentum transfers the exponential in the square bracket, 



Eq. (|A3|), simply vanish leaving a factor of unity which leads to a much reduced entropy production 
rate decreasing quadratically with time. The long range interactions (small momentum transfers) are 
therefore very important in expanding plasmas and sensitive to screening or the cutoff as is applied in 
some models ra. 



[1] J. P. Blaizot, A.H. Mueller, Nucl. Phys. B289, 847 (1987). 

[2] K. Kajantie, P. V. Landshoff and J. Lindfors, Phys. Rev. Lett. 59, 2517 (1987); K. J. Eskola, K. Kajantie 
and J. Lindfors, Nucl. Phys. B323, 37 (1989); K. J. Eskola, K. Kajantie and J. Lindfors, Phys. Lett. B 
214, 613 (1989). 

[3] K. Geiger, Nucl. Phys. A566, 257c (1994). 

[4] T. S. Biro, E. van Doom, B. Miiller, M.H. Thoma, X.-N. Wang, Phys. Rev. C 48, 1275 (1993). 
[5] K. J. Eskola, X.-N. Wang, Phys. Rev. D 48, 1275 (1993). 
[6] L. Xiong and E. Shuryak, Phys. Rev. C 49, 2207 (1994). 
[7] J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 
[8] E. Shuryak, Phys. Rev. Lett. 68, 3270 (1992). 
[9] G. Baym, Phys. Lett. 138B, 18 (1984). 
[10] S. Gavin, Nucl. Phys. B351, 561 (1991). 

[11] P. Danielewicz and M. Gyulassy, Phys. Rev. D 31, 53 (1985). 
[12] A. Hosoya and K. Kajantie, Nucl. Phys. B250, 666 (1985). 
[13] S. Gavin, Nucl. Phys. A435, 826 (1985). 

[14] G. Baym, H. Monien, C. J. Pethick, and D. G. Ravenhall, Phys. Rev. Lett. 64, 1867 (1990); Nucl. Phys. 
A525, 415c (1991). 

[15] G. Baym, H. Heiselberg, C. J. Pethick and J. Popp, Nucl. Phys. A544, 569c (1992); H. Heiselberg, in 
proceedings of Pre- equilibrium Parton Dynamics in Heavy Ion Collisions, Berkeley, 1993, LBL-34831, p. 125, 
Ed. X.-N. Wang. 

[16] H. Heiselberg, Phys. Rev. D 49, 4739 (1994). 

[17] M. H. Thoma, Phys. Rev. D 49, 459 (1994). 

[18] For relativistic particles the bulk viscosity vanishes to order o? s in the coupling constant (see Ref. |l^ , [h}[|j7| ) . 
Higher order contributions arise from the trace anomaly (R. Pisarski, private communications). 



14 



[19] M. Gyulassy and X.-N. Wans, Nu cl. Phys. B420, 583 (1994); X.-N. Wang, M. Gyulassy and M. Pliimer, 



LBL-35980, 1994, |hep-ph/9408344| , Phys. Rev. D., in press. 
[20] A. Weldon, Nucl. Phys. A566 581c (1994). 
[21] A. Weldon, Phys. Rev. D 26, 1394 (1982). 

[22] M. Gyulassy, A. Selikhov, Phys. Lett. B316, 316 (1993); H. Heiselberg, Phys. Rev. Lett. 72, 3013 (1994). 

[23] X.-N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 ( 1991). 

[24] P. Levai, B. Muller and X.-N. Wang, LBL-36594, 1995, frep-ph/9412352 . 



[25] K. J. Eskola, K. Kajantie, P. V. Ruuskanen, Nucl. Phys. B323, 37 (1989); HU-TFT-94-6 (1994). 
[26] X.-N. Wang, LBL-36885, 1995, in proceedings of 11th International Conference on Ultrarelativistic Nucleus- 
nucleus Collisions, Monterey, Jan. 9, 1995, edited by A. M. Poskanzer, J. W. Harris and L. S. Schroeder. 

FIG. 1. Total entropy in one-dimensional viscous hydrodynamical expansion of Eq. (M) for various viscosities 
n — fjT 3 . For comparison, the upper limit on the entropy by solving the Boltzmann equation, Eq. (^), is shown 
with dashed curve for 9 = t v in the case fj = 16. In both cases To = h/ro is assumed. 

FIG. 2. Qualitative picture of the expanding system following a nuclear collision. The formation time from 
— to, free streaming from tq — 9, and hydrodynamical flow for r>9. 

FIG. 3. The function g(r(x)) found by solving Eq. ([Hi]) numerically. Its dependence on the relaxation time 
9 = 9o(t/to) p is shown for various 9q and p. 

FIG. 4. Entropy production in the relaxation time approximation as function of the relaxation time 
9 = 9q{t/tq) p for various powers p. The entropy found by Baym |9J (see text) is shown by dashed curve. 
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